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ABSTRACT 

A novel experiment to study the interfacial shear instability between two liquids is described. Density- 
matched immiscible liquids are confined between concentric cylinders such that the interface is parallel to the cylinder 
walls. Interfacial waves that develop because of viscosity difference between the shearing fluids are studied as a func- 
tion of rotation rate and depth ratio using optical techniques. Conditions of neutral stability and the most unstable 
wavenumber agree reasonably well with predictions from linear stability analysis of the Navier-Stokes equations. 
Illumination using laser sheets allows precise measurement of the interface shape. Future experiments will verify 
the correctness of weakly nonlinear theories that describe energy transfer and saturation of wave growth by nonlinear 
effects. Measurements of solitary wave shapes, that occur far above neutral stability, will be compared to similar 
measurements for systems that have gravity as an important force to determine how gravity effects large distur- 
bances. These results will be used to interpret slug and annular flow data that have been obtained in other pg stud- 
ies. 

INTRODUCTION 

Cocurrent flows of immiscible fluids are found in process flows, (e.g., condensers, oil-water separators, 
iron-slag separators), hydrocarbon production and transportation pipelines and there may be future opportunities to 
replace contacting devices that rely on emulsification (e.g., liquid extractors) with separate flowing phases that have 
high interphase transport coefficients. Critical information necessary for design and operation of these devices is 
whether the phases will remain separate (perhaps due to density differences) or if mixing will occur. While 
prediction of the phase configuration has been done using flow maps [1] of dimensional variables and stability 
analysis has been done with constitutive criteria [2, 3, 4, 5], the validity of these procedures is generally confined to 
conditions close to the conditions of the experiments on which they were based. Further, the true predictive 
capabilities of these methods are lacking. If a base state with separate phases can be defined, a rigorous linear 
stability analysis that employs numerical solutions of the governing equations [6,7] can be done to determine if the 
base state is stable. While linear stability analysis can predict the conditions when waves will form and grow with 
distance, nonlinear theories are needed to determine if saturation of these waves will occur [8,9]. These rigorous 
nonlinear theories are valid only if the wave slope remains small; if the waves get too large, no predictions can be 
made. Consequently, experiments are needed to determine the long distance evolution of disturbances. 

Unfortunately, conduits available in laboratories are often too short and industrial pipelines too irregular to determine 
definitively what the long distance state of a linearly unstable flow is. 

To solve this problem, we have devised a new experiment that uses a two-layer flow of density matched flu- 
ids confined between concentric cylinders. The interface is parallel to the cylinder walls and rotation of the outer 
cylinder, for fixed inner cylinder, induces shearing that can lead to instability. The waves that form are completely 
analogous to interface waves formed by a shear flow in a flat channel. For the rotating device these have crests that 
are parallel to the rotation axis and exhibit interface deformation in the radial direction (r) that is initially periodic in 
the direction of rotation (0). This device can simulate extremely long distance evolution of disturbances by simply 
waiting a sufficiently long time. For example, experiments shown below indicate that waves with significant 
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amplitudes will form only after several minutes at rotation rates large enough to cause instability. This corresponds 
to over 100 meters of length — a distance not usually available for a well-controlled experiment. 

Interfacial wave growth curves for systems such as gravity driven films, pressure driven pipe or channel 
flow and two-layer Couette flows have two distinct shapes. Figure 1 shows temporal wave growth curves calculated 
numerically [10] using a two-layer Orr-Sommerfeld analysis and full linearized boundary conditions. In figure la the 
region of positive growth is bounded away from 0 wavenumber. Figure lb shows the second shape where the 
growth rate has a maximum away from 0 but is positive all the way down to 0 once any wavenumber becomes un- 
stable. The primary question suggested is: what generic differences are there in the very long time evolution of 
waves for the two different types of growth curves? Previous work [11] suggests that wave modes can be generated 
in linearly stable regions by nonlinear interactions, but it is not clear if they stay at small amplitudes or grow into 
the dominant waves. Gravity driven films, once unstable, always look like figure lb, channel flows can exhibit ei- 
ther of the two behaviors. 

The two-layer Couette geometry changes from la to lb as h, the ratio of the inner (less viscous film) layer 
thickness/gap thickness, is altered. Thus it is possible to compare the behavior of both types of growth curves 
using well-controlled long time experiments. In this paper we describe the device, the experimental procedures and 
give some preliminary results on regions of instability and wavelengths of observed waves. 


EXPERIMENTAL 

f 

Wave measurements have been made in a rotating concentric cylinder test cell. The diameter of the inner 
cylinder is 19.50 cm and the diameter of the outer cylinder is 21.50 cm resulting in a gap thickness of 1.00 cm. 

The gap width is 5 cm. The cylinders are loaded into a Weissenberg rheogoniometer such that the inner cylinder is 
fixed, the outer cylinder may be rotated, and the axis is aligned with gravity. The concentricity of the cylinders is 
aligned to a tolerance of about 75 fvn while the tilt of the cylinders is aligned to a tolerance of about 25 /m. 

To perform measurements in our system, the outer wall is transparent to permit direct imaging of the inter- 
face. We use several lighting techniques to probe the experiment — white light, a vertical plane of laser light, and a 
horizontal laser sheet. The two laser methods provide useful quantitative data that direct visualization with white 
light does not allow. Figure 2a is a diagram of the vertical laser setup. It shows a plane of laser light projected on 
to the interface at a large incident angle. The laser light is visible at the outer surface of the outer wall and at the in- 
terface. The distance between the two laser images is proportional to the fluid depth ratio. Further, the image on the 
interface is a profile of the interface shape. Figure 2b is a diagram of the horizontal laser setup. The horizontal laser 
set up is designed to show a wave tracing at a constant vertical position. 

The two fluids used for the experiments whose results are presented here are Dow 710, which is a 
phenylmethyl polysiloxane fluid and a mixture of ethylene glycol, water, and Pink Bismuth - an Osco® brand upset 
stomach remedy. The Pink Bismuth is used as a source of refractive particles. The viscosity and density of Dow 
710 are 0.555 Ns/m? and 1110 kg/m^ respectively. The viscosity and density of the ethylene glycol solution are 
0.0151 Ns/m 2 and 1 108 kg/m? respectively. Dow 710 is loaded on the outer cylinder after the ethylene glycol 
mixture has been added to the channel. 

To map the stability boundary of our system, the velocity of the outer plate was set to about 0.05 m/s and 
periodically increased until waves formed. The exact value of the initial velocity is dependent upon the stability of 
the interface. If conditions change much from the conditions when the fluids were density matched, larger velocities 
are needed to maintain the interface. 
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RESULTS AND OBSERVATIONS 


The growth curves in figure 1 predict that waves will form if the rotation velocity is sufficiently large and 
become detectable after a sufficiently long time. Figure 3 is a plot of Reynolds number versus dimensionless inner 
film thickness. The three curves correspond to dimensionless growth rates of 0.0, 0.001, 0.015. The circles show 
the experimentally - observed onset of waves. For a height of 0.7, waves were observed very close to the predicted 
onset conditions. A growth rate of 0.001 corresponds to a time constant of approximately 5 minutes and a flow dis- 
tance of 100 m! For smaller heights, a larger growth rate was required to observe waves. This may have been be- 
cause we did not wait long enough for waves to form, or it could be because nonlinear effects caused saturation at 
amplitudes too small to observe. Further experiments and nonlinear analysis are needed to resolve this issue. 

Observed wavenumbers match well with the predicted fastest growing waves. Figures 4a-d are plots of the 
theoretical growth rate curves for the most unstable wave mode and experimental data. Wavenumbers for experi- 
mental data are depicted with a dashed line. Most of the stability points lie near the maximum growing wave num- 
ber. The wavenumber for figure 4b lies closer to 1/2 the most unstable wavenumber. This could have been excited 
when the rotation rate was lower and thus the most unstable wave mode occurred at a lower wavenumber. This 
behavior has been seen in previous experiments with conditions such that the system was unstable for all Reynolds 
numbers. In those experiments, a small amplitude long wave was excited at a low plate velocity and grew in 
amplitude with increasing velocity. An alternative explanation is that there could have been nonlinear interactions 
that caused transfer of energy to lower wavenumber. Further experiments that carefully watch the time evolution of 
the wavenumber will be done to see which scenario is occurring. 

CONCLUSIONS 

Two-layer shear flow interfacial instability can be studied in the geometry of a rotating Couette flow that al- 
lows observation of wave evolution at times and (effectively) distances much longer with much better control than is 
possible in channel flows. Density matching of the fluids, necessary to create the vertical interface geometry, essen- 
tially removes gravity as an important force. By doing long time experiments, it is possible to examine the impor- 
tance of modes near zero wavenumber that have much lower growth rates than shorter modes. These experiments 
will enable careful checks of nonlinear theories that predict energy transfer between modes and saturation of ampli- 
tudes at sufficiently long times. Direct white light imaging using a standard video camera allows measurement of 
the wavenumber and wave speed. The precise interfacial wave shapes, which are particularly important away from 
neutral stability where the shape is nonsinuous, can be accomplished using laser sheets to illuminate the flow cell. 

It is hoped that these measurements of nonsinuous waves will help to interpret annular and slug flow behavior that 
is observed in (xg studies that have been done. 
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figure 2a: Vertical laser setup. 
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figure 2b: Horizontal laser setup. 
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figure 3: Stability curves for three characteristic growth rates plotted with experimental data 
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